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ABSTRACT 


It  is  shown  how  a  reliability  shorthand  can  be  imple¬ 
mented  on  a  handheld  calculator. 

Assuming  constant  failure  rates,  basic  structures  are 
used  to  show  how  the  shorthand  can  be  applied.  Several  exam¬ 
ples  are  worked  out  that  show,  how,  with  component  failure 
rates  as  input,  a  handheld  calculator  can  be  used  to  compute 
the  reliability  of  a  system. 

Two  TI-59  programs  are  provided  as  a  computational  aid. 
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I.  INTRODUCTION 

Systems  and  coaponents  can  be  in  either  of  two  states: 
either  they  are  functioning  or  they  have  failed.  The  abil¬ 
ity,  that  a  system  stays  functioning  over  a  predetermined 
time  interval  is  called  its  reliability.  It  is  generally  not 
realistic  to  assume  that  a  system,  say  a  lightbulb,  will 
fail  at  a  specified  time,  but  rather  that  T,  the  time  to 
failure,  is  a  random  variable  which  has  a  probability  dis¬ 
tribution  that  can  be  specified.  The  probability  distribu¬ 
tion  for  a  time  to  failure  is  called  its  life  distribution. 
In  this  paper  we  will  solely  be  concerned  with  one  specific 
type  of  life  distribution  which  is  especially  important  in 
reliability  theory  and  practice,  the  exponential  distribu¬ 
tion.  It  has  the  property  that  the  remaining  life  of  a  used 
component  is  independent  of  its  age  (the  "memory less”  prop¬ 
erty)  ,  i.e.  a  functioning  component  is  always  as  good  as 
new,  the  failure  rate  is  constant.  The  memoryless  property 
is  the  basis  for  a  reliability  shorthand,  one  that  can  be 
implemented  on  a  handheld  calculator. 

Depending  on  the  size,  structure  and  life  distribution 
of  a  system,  probability  statements  about  its  time  to 


8 


failure  are  in  general  not  easily  achieved.  Porming  the  sua 
of  independent  life  lengths  (i. e.  convolving  the  corres¬ 
ponding  life  distributions!  requires  knowledge  of  integral 
calculus  and  coaputations  can  becoae  rather  tedious. 

In  zhe  case  of  the  exponential  distribution,  though, 
coaputations  can  be  siaplified  by  translating  the  problem 
into  a  simple  shorthand  notation  and  using  this  shorthand  as 
input  for  soae  computing  device. 

In  this  paper  we  will  show  how  a  reliability  shorthand 
can  be  implemented  on  a  handheld  calculator.  Basic  struc¬ 
tures  are  used  to  show  how  the  shorthand  can  be  applied.  Two 
TI-59  programs  are  provided  as  a  computational  aid.  Formulas 
for  the  convolution  of  up  to  four  exponential  random  varia¬ 
bles  can  be  found  in  Appendix  A.  Appendix  B  contains  a  user 
guide  to*the  TI-59  programs. 
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II.  JHE  CONCEPT  2£  h  H ELIABILITY  SHORTHAND 

A.  BASIC  NOTATION 

The  survival  function  of  a  life  length  can  be  derived 
from  the  distribution  function. 

Let 

T  :  life  length 

?(t)  =  P(T<t)  be  the  distribution  function  of 

Then 

?  (t)  ■  P  (T>t) 

=  1-F(t) 

is  the  survival  function  of  T. 

-It 

In  the  case  of  the  exponential  distribution,  F(t)=  e  , 
where  h  is  the  failure  rate.  Translated  into  shorthand,  the 
life  distribution  is  denoted 
EXP  (  K  ) . 

B.  CONVOLOTION  OF  DISTRIBOTIONS 

Bher.  independent  random  lives  are  summed  up,  the  corres¬ 
ponding  life  distributions  have  to  be  convolved  to  determine 
the  probability  that  the  sum  of  the  lives  will  exceed  a  spe¬ 
cified  time  t.  Let 

T,  ,T,  :  independent  life  lengths 
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F*  (t)  ,Tt  (t)  :  the  correspo na ing  survival 

f u  nctions 

f-t  (t)  ,ft  (t)  :  the  corresponding  density  functions 
T  =  T.,  +  T*  :  the  total  life  length 

Then 

F  (t)  =  P  (T>t) 

=  P  (T^  ♦  T z  >t) 
t 

■  F,  (t)  +  j  7*  (t-s)  f<  (3)  ds. 
o 

This  means  that  T  will  exceed  a  specified  *-ime  t  when 
-either  T^  exceeds  t 

-or  T,  is  smaller  than  t,  say  equal  to  s,  and  Ta 
exceeds  t-s. 


Integration  with  respect  to  s  (i.e.  summing  over  all  possi¬ 
ble  values  of  s  )  is  called  the  convolution  of  T,  and  T2 . 
When  T«  and  Tz  are  both  exponentially  distributed  with  fai¬ 
lure  rates  and  ,  i.e. 

F«  (t)  =  e”*** 

Fj  (t)  =  Q1**  , 
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This  shorthand  notation  is  heuristic* lly  apparent.  We  can 
visualize  a  1  component  /  1  spare  system  with  EXp(/t,4)  and 
Exp  ( /i,  )  lives  respectively.  From  component  1  the  system 
has  an  EXP  ( )  life  to  begin  with.  When  component  1  fails, 
the  system  has  an  extra  EXP  (  Az)  life. 

C.  MIXTURE  OF  DISTRIBUTIONS 
1  •  MIX-Notation 

In  the  previous  chapter,  we  formed  the  sum  of  inde¬ 
pendent  random  lives,  which  each  had  weight  one,  i.s. 

T  =  +  T2  . 

Now  consider 

T^  with  probability  p^ 

T  =  < 

^  Tj  with  probability  p£ 
where  p.,  +  P*  =  1  . 

Let  D.,  and  Dj,  be  the  probability  distributions  of  the  random 
variables  T^and  Tj  respectively.  The  corresponding  survival 
functions  are  P,,  (t)  and  F2  (t)  . 

Then 

F(t)  =  p,  P,  (t)  +  Fa  (t)  . 

In  shorthand,  the  mixture  of  distributions  and  Dz  with 
respect  to  the  mixing  probabilities  p0  and  p2  is  denoted 
MIX  [  p^  ,  p2  Oj  ]. 
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2 •  Distributive  Law 
How  let 


T  =  Tj  ♦  T* 


where 


r 


T1  with  probability  p 


T '  -  \ 


Tt  with  probability  1-p. 


Then 


r 


T  =  T3  ♦  \ 


T4  with  probability  p 


r 


^  T2  with  probability  1-p. 


Tj  +  T.,  with  probability  p 


T  =  < 


^  Tj  +  T2  with  probability  1-p. 

The  distributive  law  holds  due  to  the  fact  that  the  sum  of 
the  mixing  probabilities  for  T.,  and  r2  is  equal  to  one. 

The  survival  function  of  T  can  be  found  by  convolution: 

* 

P(t)  =  F,  (t)  +  J(  p?4  (t-s)  *  (l-p)fj  (t-s)  )  f,  (s)ds. 
o 

With  D„  ,  D2  ,  Dj  being  the  probability  distributions  for  T1r 
Tj ,  Tj  t  the  distributive  law  can  be  applied  to  the  shorthand 
notation: 

D  j  ♦  HIX  [pD1  ,  (1-p)  D2]  =  HIX  [  p((D„  ♦  Dj )  , 

(i-p)  p2  +  D?)  ]. 
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Graphically  this  can  be  represented  as  follows 


Figure  1:  Distributive  Property  of  the  MIX-Notatior. 

3-  Degen eracy  at  the  Origin 
Let 

P  (T  =0 )  *  1. 

Then  the  distribution  of  r  is  degenerate  at  zero. 

In  shorthand  notation,  such  a  distribution  is  called  the 
ZERO-distr  ibution. 

Now  let  T  =  +  T0 

where  T*  and  T„  have  probability  distributions  D,,  and  ZERO 
and  survival  functions  f^(t)  and  ?0  (t)  respectively. 

Then 

t 

F{t)  =  Pi  (t)  +  )F0  (t-s)  %  (s)  ds 
o 

=  F„  (t )  . 


i 


The  ZEBO-distribution  doesn't  add  anything  to  another  dis 
tribution,  so  for  instance 


D4  ♦  ZERO  = 

D2  +  HIX[  pD,,  ,  (1-p)  ZERO  ]  =  HIX[p(Dyf+  D2),(1-p)D2 
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III.  APPLYING  A  RELIABILITY  SHORTHAND 
After  this  brief  survey  over  the  concept  of  a  reliabil¬ 
ity  shorthand  we  will  now  show  how  the  shorthand  can  be 
applied.  To  do  sc  we  will  use  basic  structures.  Part  A  of 
this  chapter  will  give  examples  whose  representation  in 
shorthand  requires  only  basic  notation  described  in  Chapter 
II,  Parts  A  and  B,  whereas  Part  B  of  this  chapter  will  give 
examples  whose  representation  in  shorthand  makes  use  of  the 
MIX-notation  and  the  ZEHO-d istr ibution . 

A.  SUMS  OF  EXPONENTIALS  WITH  WEIGHT  ONE 
1 •  Simple  Series  Syst e m 

A  series  system  is  a  system  which  is  functioning, 
when  all  its  components  are  functioning.  A  two-component 
series  system  can  be  graphically  represented  as  shown  in 
Fig .  2  . 

Let 

T  :  life  of  the  system 

Tf  :  life  of  component  1 

Ta :  life  cf  component  2 

F4  (t)  =  survival  function  of  component  1 
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in 

n 

D 

Pigure  2:  Two-Component  Series  System 


Then 


F*  (t) 


survival 


function  of  component  2 


T  =  min  (  T., ,  T2) 

F  (t)  =  survival  function  of  the  system 
=  P(  min  (T4,  Tj)  >  t) 

-  P(  T^>  t,  T2>  t) 

Assuming  independence  of  the  two  components 
f (t )  -  p (  t^ >  t)  P(  r2>  t) 

=  (t)  ?z  (t) 

-a,  £ 

=  e 

The  shorthand  notation  for  this  system  is 
EXP  (/i,,  +  /t2)  . 
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This  is  intuitively  apparent,  as  the  system  has  an  exponen¬ 
tial  survival  function  with  failure  rate  . 

2.  Simple  Parallel  System 

A  parallel  system  is  a  system  which  is  functioning, 
when  at  least  one  of  its  components  is  functioning.  A  two- 
component  parallel  system  can  be  graphically  represented  as 
follows: 


Figure  3:  Two-Component  Parallel  System 
Let 

Tn  ~  EXP  (  i  )  ,  T1  ~  EXP  {  l  )  . 

Then 

T  =  max {T^ ,  T*) 

F(t)  =  P  (ma  x  (T, ,  Tt)  >  t) 

=  1-  P  (max  (T,,  ,  r, )  <  t ) 

*  1  -  P(T,<  t,  t) 
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Assuming  independence  of  the  two  components, 
F(t)  =  1  -  Ptt,  *  t)  P  (T,<  t) 


*  1  *  <t>  f2  (t) 


.  .j  “(it  .  ..  . 

*  1  -  (1-e  )  (1-e  ) 


,  ,,  -At  -2At  % 

=  1  -  (1-2?  ♦  e  ) 


-  -Art  -7 At 

-  2e  -  e 


The  shorthand  notation  for  the  system  is 
EXP  (2  K  )  *  EXP  (A,)  . 

This  follows  intuition  as  the  system  has  an  EXP(2/l)  life  to 
begin  with  and  when  one  component  fails  it  has  an  extra 
EXP  (  i  )  life  due  to  the  mamoryless  property  of  the  exponen¬ 
tial  distribution. 


3.  Standby-system  with  Dissimilar  Components 

Suppose  a  system  consists  of  two  components,  one 
active  and  one  spare.  The  active  component  stays  in  service 
until  it  fails  and  then  immediately  is  replaced  by  th® 


spare. 

Let  the  time  to  failure  of  the  two  components  be 
T,  ~  EXP  (/i,,  )  and  Ta  ~  EXP  (  )  respectively. 

Then  the  system  time  to  failure  is 
T  »  T<,  ♦  Ta 

and  the  survival  function  of  the  system  is 


F  (t )  *  P  (T  >  t ) 


The  shorthand  notation  for  the  system's  survival  function 
should  be  obvious.  The  system  has  an  EXP(  /i4 )  life  from  the 
active  component  and  an  additional  EXP(/ij)  life  from  the 
spare.  Sc  the  shorthand  notation  is 
EXP  {  )  +  EXP  (  )  . 
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B.  SUMS  OF  EXPONENTIALS  WITH  WEIGHT  BETWEEN  ZERO  AND  ONE 


The  examples  given  in  the  previous  chapter  only  involved 
exponential  lives  with  weight  one.  Now  we  will  look  at  some 
structures,  whose  survival  function  has  a  shorthand  notation 
which  includes  the  MIX-notation  and/or  the 
ZER  O-distr ibution . 

1 .  Parallel  System  wit h  Pis si nil ar  Failure  Rates 

The  notion  of  a  parallel  system  has  been  introduced 
in  Chapter  III. A. 2  .  We  now  look  at  the  case  where 
T^EXPti,,)  and  T4  EXP  (  /tj  )  . 

Than 


T  =  max  (Tn ,  T,) 

F(t)  »  P(  max  (T^ ,  T2  J  >  t) 

«  1-P(  max(T^  ,  T,)  *  t) 
=  1-P(T^  t,  T,$  t) 


Assuming  independence  of  the  two  components 
F  ft)  »  1-P(T,$  t)  P(T2<  t) 

*  1-fc,  (t)  ?t  (t) 

=  (i-e"^  -e" } 

-^4 t  -  j  t  ■*  l  A*  *  )£ 

=  e  ♦  e  -a  . 
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To  find  the  shorthand  notation  of  the  system  consider  all 
the  ways  which  lead  to  the  survival  of  the  system: 

-either  both  components  survive 

-or  component  1  fails  and  component  2  survives 

-or  component  2  fails  and  component  1  survives. 


If  one  component  fails  and  one  survives,  in 


fraction 


of  the  cases  the  survivor  will  be  component  1  and  in 
fraction  of  the  cases  it  will  be  component  2. 

This  can  graphically  be  represented  as 


EXP  ( ) 


A'a+A'  i 


EXPtX,,  +  it)  ♦ 


t  4 1 


EXP  (  i ,  ) 


Making  use  of  the  Mix-notation  the  shorthand  notation  then 


EXP  Ui  +  / iz)  *  MI  X[  EXP  (  EXP  (  )  ] 

and  using  the  distributive  property  it  becomes 

MIX[  il—tEXP  ( K^)  *  EXP(  4  +  /tj)  )  , 

- (EXP(/tj)  ♦  EXP(  +  It)  )  ]. 

Am  *  A>x 
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As  a  check  to  see  that  this  shorthand  notation 


represents  the  survival  function  of  the  system,  we  derive 


the  survival  function  from  the  shorthand  notation: 


F(t)  = 


Ax 


-A*t  f 

A*  C 


1 * * j 


’  ds} 


-Ait  ~(A/s+At)t 

=  e  +e  -  e 

This  verifies  that  the  shorthand  notation  indeed 
represents  the  system’s  survival  function. 

2 .  Series  System  with  One  Spar e 

Let  us  now  look  at  a  two-component  series  system, 
whose  components  have  dissimilar  failure  rates  with  one  com 
ponent  having  a  spare: 

Component  1  has  the  constant  failure  rats  and 
component  2  and  the  spare  have  the  constant  failure  rate  Lz 
The  spare  can  only  replace  component  2. 

Let 

F,,  (t)  :  the  survival  function  of  component  1 
Pj  (t)  :  the  survival  function  of  the  standby 
system  component  2  with  its  spare. 


23 


24 


A, 


■fraction  of  the  tine  if  will 


If  one  component  fails,  in 

A, 

be  component  1,  which  means  that  the  system  will  not  sur- 

iz 

vive;  in - - fraction  of  the  time  the  failing  component 

A*+A 1 

will  be  component  2. 


This  can  graphically  be  represented  is 


ZERO 


EXP  (  l,  +  )  ♦ 


exp  (  A*  +  Az ) 

Using  the  MIX-notation  the  survival  function  then  is 

EXP  (  ♦  iz  )  +  Mixr  ^7-  ZERO,  -;— ■*<— EXP  (  iA  +  Az  )  ] 

=  MIX[  ■  -(ZERO  +  EXP (  A*  +  Az  ))  , 

A*  Az 

— (exp  (i,  +  / 12 )  +  EX?  (  A<  ♦  At )  )  ] 

A*  *  Az 

-  MIX[-_^T-(EXP(  a,  +  At))  , 

A*  +  At 

— ^k_-(EXP  (An+ii)  +  EXP(  h*lt)  )  ]. 

ziit  +  Al 

To  prove,  that  the  shorthand  notation  does  represent  the 

survival  function,  we  derive  the  latter  from  the  shorthand: 

-(A*  +  Ax)t  .  -(^  +  4,)* 

1 — t“  £  +  *7 — — —  [  e  + 

A  a  *  A  ^  A*\  +  At 


P  (t) 


J  -  ( _s)  -  ( A,*  +  A.\ 

J  t  U-.  +  'U)  t 


.  -[A«  +  k \)t 

=  (  1  +  /t2t  )  e 


This  is  the  previously  found  result  and  this  verifies,  that 
the  shorthand  notation  does  represent  the  system's  survival 
function. 

3.  Two-out-of-Th ree  S  y  stem 

As  a  last  example  in  this  chapter,  we  will  look  at  a 
Two-out-of-  Three  system. 

Consider  a  three  component  system,  whose  components 
have  constant  failure  rates  /i4  ,  /t2and  /tj  respectively.  The 
system  is  functioning,  as  long  as  two  out  cf  three  compo¬ 
nents  are  functioning  (see  Fig.  6)  . 

In  other  words,  the  system  is  functioning  as  long  as 
there  is  a  path  through  the  system  . 

Alternatively,  the  system  can  be  visualized  as  a 
par allel -series  system  (compare  Fig.  7). 

The  survival  function  of  the  system  is 

F(t)  =  P  (Ty,  >t  A  >t)  +  ?  (T,  >t  A  Ts  >t) 

+  P (Tj  >t  A  Tj  >t) 


P(  (T, 

>t 

A 

>t) 

A 

>-  A 

** 

>t)  ) 

P(  (T^ 

>t 

A 

>t) 

A 

(Ta 

>t  A 

>*-)) 

P(  (T< 

>t 

A 

*3 

>t) 

A 

(t2 

>t  A 

T3 

>t)  ) 
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Thus 

?<t)  =  P  (T^  >t  a  Ta  >*)  ?  (T.,  >t  a  T,  >t) 

♦  P(T3  >t  A  T3  >t) 

-  3P  (T^  >t  A  Tt  >t  A  T3  >t) 

♦  P  (T^  >t  A  T2  >t  A  r3  >t) 

Therefore,  and  assuming  independence  of  rhe  components, 

?(t)  =  P{T^  >t)  P  (T2  >t)  +  P  (T^  >T)  P(T3  >t) 

♦  P(Ta  >t)  P(T3  >t) 

-  3P  (T^  >t)  P(T2  >t)  P  (Tj  >f) 

♦  P  (T^  >t)  P(T2  >z)  P  (Tj  >t) 

=  P  ( T,,  >t)  P  (T2  >t)  +  P(T^  >t)  P(T3  >t) 

♦  P(T2  >t)  P  (Tj  >t) 

-  2P  (T^  >t)  P(T2  >t>  P(T3  >t) 

-  (A*  +/tj  )  f  *  -( 

—  3  ▼  e  r  e 

~  (  4,*>  +  Ai  +  4,i )  t? 

-  2e 

Now  let  us  consider  all  the  possible  ways,  in  which  the  sys¬ 
tem  can  survive: 

-  either  all  components  survive 

-  or  component  1  fails  and  component  2  and  3 
survive 

-  or  component  2  fails  and  component  1  and  3 
survive 
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-  or  component  3  fails  and  component  1  and  2 


survive. 

If  a  component  fails  and  the  other  two  survive,  in  — — — - 

XA  +  A,  +A 

fraction  of  the  time  it  will  be  component  i  ,  i  =  1,2,3  . 
This  can  graphically  be  represented  as 


EXP  (A,  ♦  A*  * 


The  shorthand  notation  then  is 


EXP  (  A*  +  Ax  ♦  AS)  *  MIX[ 


- - ^ — -—EXP  (  Az  +  A,  )  , 

A*  +  Ai+ Aj 

- - ^ - —EXP  (  A„  ♦  A,  )  , 

A„  +  Kx  tA3  3 


As 


A*  +■  +A j 


EXP (  A„  +  Av)  ], 


=  HIX( 


t  Ax  +  A3 

A% 

A*  +  Ax  +  Aj 

_ iLi _ 

At  *  A%  +A3 


(exp(  Ax  *  A3 ) 
(ex p  (  Aa  +  A,) 
(ex p (  A„  *  Ax) 


♦  exp(A„*  Az  +  A3  )) , 

+  exp  (  Aa  *  Ax  *  Aj  ) ) , 

♦  exp  ( At  *  Az  +  A3 ) )  ]. 
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Again,  as  a  check  that  the  shorthand  notation  represents  the 
survival  function,  let  us  derive  the  survival  function  from 
the  shorthand  notation: 


K  a  r 

F(t)  =  m — j —  t  * 

A**  A » +  X  i  ± 

+  j  e  M*+/ l»)<? 


Xz  r  -  (  X«,+X*)  £ 

♦  -  t 

X„  +  x2  +  x3  l  , 

I  - (X-wX*-# X*) £4-5)  -(X<|tXj)s 

4  J*  U„U,)£  dsj 


X*  +  Xi  +  x? 


(  -C>i4+'lz+Xi)  -(Aa+Az^S 

+  )C  (A<+A  i)e  •  eU] 


-(Ai  +  As  )£  -(A*  +  A-$)i  -  ( X„  +  X*)  £ 

=  £  +  £  *  t 


-it 


-  (  X*  +  X*.  +  X  j )  £ 


The  result  again  proves  that  the  shorthand  notation 


indeed  represents  the  survival  function  of  the  system. 


IV.  Ij1£LEMENTING  T  HE  SHORTHAND  ON  THE  TI-59 
The  concept  of  a  reliability  shorthand  is  introduced  in  ^he 
course  "Reliability  and  Weapons  system  Effectiveness 
Measurements",  OA  4302,  at  the  Naval  Postgraduate  School, 
Monterey.  Most  students  talcing  the  course  are  in  the 
Operations  Research  (OR)  -  Curriculum. 

The  choice  of  the  TI-59  as  the  computing  device,  on 
which  the  shorthand  was  to  be  implemented,  was  based  on  the 
fact,  that  each  student  in  the  OR-Curri culum  is  issued  a 
TI-59  for  use  in  basic  probability  and  statistics  courses. 
Thus  almost  every  student  at  the  Naval  Postgraduate  School, 
who  is  introduced  to  the  shorthand,  is  familiar  with  the 
TI-59  and  has  access  to  such  a  calculator. 

A  program,  that  uses  the  shorthand  notation,  times  to 
failure  and  failure  rates  as  input,  should 

-  calculate  the  survival  probability  of  basic 
structures  /  small  systems  and 

-  require  moderate  computation  time. 

To  achieve  these  requirements  it  was  decided  to  incorporate 
all  solutions  for  the  convolution  of  up  to  four  exponential 
random  variables  in  the  program.  The  formulas  that  were  used 
are  given  in  Appendix  A. 
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Two  programs  are  provided  in  this  paper. 

Program  1  can  be  used  when  all  rates  are  dissimilar  cr 
all  are  the  same.  It  uses  the  formulas  on  pages  37  and  38 
onl  y. 

Program  2  can  be  used  for  the  general  case.  It  makes  use 
of  all  the  formulas  given  in  Appendix  A.  The  program 
includes  a  sorting  routine  that  determines  the  applicable 
formula  from  the  entered  failure  rates. 

A  user  guide  to  the  two  programs  is  provided  in  Appendix 
B. 
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7. 


SUMMARY 


There  is  a  reliability  shorthand  that  denotes  the  survival 
function  of  a  system,  assuming  that  the  failure  rates  of  all 
components  are  constant. 

This  shorthand  can  be  implemented  on  the  TI-59  handheld 
calculator.  With  failure  rates,  time  to  failure  and  short- 
hem  d  as  input  the  TI-59  calculates  the  survival  probability 
of  the  system. 

Knowledge  of  calculus  is  not  necessary  to  use  this 
method,  whereas  the  standard  procedure,  finding  the  survival 
probability  by  convolution,  requires  knowledge  of  integral 
calculus. 

The  choice  of  the  TI-59  as  the  computing  device  for  the 
implementation  of  the  shorthand,  though,  implied  limita¬ 
tions;  the  number  of  failure  rates  is  limited  due  to  the 
limited  storage  capacity  of  the  TI-59,  and  computing  times 
are  comparatively  long.  The  TI-59  can  therefore  only  be  used 
for  smaller  systems,  preferably  for  the  solution  of  class¬ 
room  problems. 

For  the  solution  of  larger  problems,  the  shorthand 
should  be  implemented  on  a  state-of-the-art  personal 


33 


computer  using  a  general  algorithm  for  the  convolution  of 
any  number  of  exponential  random  variables. 
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APPENDIX  A 


CONVOLUTION  FORMULAS 

Appendix  A  contains  formulas  for  the  convolution  of  up 
to  four  exponential  random  variables. 

For  the  two  special  cases,  when  all  random  variables 
have  the  same  failure  rate  and  all  have  different  failure 
rates,  general  formulas  for  the  convolution  of  any  number  of 
exponential  random  variables  are  given. 

These  formulas  are  used  in  the  two  TI-59  programs  pro¬ 
vided  in  Appendix  B. 
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■"  .  . . .  .. 


System: 


Shorthand:  2XP(/l) 


Survival  Function:  F(t) 


-Jit 


Description : 

A  single  active  component  with  constant  failure  rate  /{. 
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System: 


Shorthand:  EXP(  h  )  ♦  EXP  (  i  )  +  ...  +  EXP  (  Ji  ) 


Survival  Function:  F(t) 


ut)°  at)*  a*)w 

— — + - + ...  + - 

0!  'll  U-4)/ 


-It 

e 


» 


i*A 


Ui)i"'  -At 
U-'IV.  e 


Description : 

A  single  active  component  with  constant  failure  rate  is 
supported  by  n- 7  identical  spares. 
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Sys  tea 


\ 


Shorthand:  EX'P  (A^)  +  E 


Survival  Function:  F(t)  = 


Is  * 


Description: 

A  single  active  component  w 
upported  by  n- 1  spares.  The  ac 
have  all  constant,  but  dissimi 
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System: 


Shorthand:  EXP(/{,„)  +  EXP (4,)  ♦  EXP(/{,a) 


Survival  Function:  F(t) 


where  A 


(B  ♦  C t) 


~  t 
e 


( /t  *  -  A,*  ) 1 


B  =  1  -  A 


Description : 

A  single  active  component  with  constant  failure  rate  /i 
is  supported  by  two  spares  with  identical  constant  failure 


System 


Shorthand:  EXP(X,4)  +SXP{/{1)  ♦  EXP  (  X.,  )  ♦  EXP  (  /L  ) 


Survival  Function:  ?(t)  = 


-l«t  2  -Ati 

Ae  +  (B  ♦  Ct  ♦  Dt  )  e 

where  A  =  - — — - = 

B  =  1  -  A 


c  =  <4,  - 


lx 


(  /iy,  “  X-l  ) 


D  = 


A, a  'Ll 


2  (  -  A*x ) 


Description : 

A  single  active  component  with  constant  failure  rate  /I* 
is  supported  by  three  spares  with  identical  constant  failure 
rate  /i2. 
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System: 


Short  hard:  EXP  (/{,„)  ♦  EXP  (/{,„)  +  EXP(/tz)  +  EXP  (  JL^  ) 


Survival  Function:  P(t) 


-  t  -A.it 

{  A  ♦  Bt) e  +  (C  +  D*) e 

—  3  A.* 


whera  A  = 


B  = 


(  'Li-A,')3 

iil 


C  =  1  -  A 
2 

A.. 

D  = 


(  Arl  ~  A.*)*’ 

-  A 

i}  lz 


U„-A.i)z 


Description : 

A  single  active  component  with  constant  failure  rate  'i* 
is  supported  by  one  identical  spare  and  tvc  spares  with  dis¬ 
similar,  constant  failure  rate  . 


System: 


Shorthand: 


EXP(^)  +  EXP(/t2)  ♦  BXP(4j)  ♦  EXP(>t,» 


Survival  Function:  F(t)  =  +  Be  ^  +  (C  +  Dt)  e 


where  A  = 


/tz  ''t'l 


B  = 


A±Al 


c  = 


x»<1  +  A, l  /ia  I  4 


D  = 


u.-Ainu^)1  u7- 

/i<l  <il 


(  A«-A  j)  (-{*  -^,3) 


Description: 

A  single  active  component  with  constant  failure  rate  i 
has  three  spares.  One  spare  has  constant  failure  rate  <ii  , 
two  spares  are  identical  with  constant  failure  rate  . 
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APPENDIX  B 


USER  GUI5S  1°  TI-59  PROGRAMS 

Appendix  B  contains  a  user  guide  to  two  TI-59  programs, 
which  use  reliability  shorthand  and  failure  rates  as  input 
to  compute  the  survival  probability  of  a  system. 

PROGRAM  1  is  designed  for  the  two  special  cases  where 
the  reliability  shorthand  is  of  the  form 

EXP  (  i  )  ♦  EXP  (  i  )  +  ...  +  EXP  (  i  ) 
or 

EXP(^)  *  EXP  (A>i)  +  ...  +  EXP  (  /ij, )  . 

In  the  first  case  the  number  of  terms  is  not  limited,  wher¬ 
eas  in  the  second  case  the  number  of  terms  is  limited  to  40 
due  to  limited  storage  capacity  of  the  TI-59.  In  this  case 
the  number  of  terms  can  be  increased  to  70  by  entering  9  in 
the  display  and  pressing  2nd  Op  1  7  . 

PROGRAM  2  is  designed  to  solve  problems  of  the  kind, 
that  were  introduced  in  chapter  III.B.  .  Due  to  limited 
memory  of  the  TI-59  the  number  of  exponential  terms  under 
one  weight  in  shorthand  notation  is  limited  to  four. 
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PROGRAM  1  :  Procedure 

1.  Use  any  library  module.  Read  in  program  1  (side  1  of 
the  magnetic  card) 

2.  Press  2nd  C'  to  initialize. 

3.  Enter  n,  the  number  of  exponential  terms  to  be  convolved, 
in  the  display  and  press  A. 

4.  Enter  time  t  and  press  B. 

5.  Enter  /i;  and  press  C  .  When  all  failure  rates  are  the 
same,  enter  -i  only  once. 

6.  a)  To  find  the  survival  probability  of  the  system, 

when  all  failure  rates  are  the  same,  press  2nd  A', 
b)  To  find  the  survival  probability  of  the  system, 

when  all  failure  rates  are  dissimilar,  press  2nd  B' . 


PROGRAM  1  :  Sample  Problems 

1.  Find  the  survival  probability  of  a  parallel  system 


( 

compare 

Chapter 

III.  A 

•  2  ) 

a) 

9 

t  = 

7  , 

n  =  2 

b) 

Short ha nd 

not  at 

ion : 

EXP  ( 

.6) 

♦  EX 

P  (-  3) 

c) 

Enter 

Com  ment 

Pr  ass 

Display 

Initi 

a  lire 

Z  1 

D 

2 

n 

A 

0 

7 

t 

B 

7 

.6 

2l 

w 

.3 

.3 

l 

r* 

.  3 

F(t) 

B' 

.2299172797 

calculation  takes  13  seconds 
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2.  Find  the  survival  probability  of  a  standby-system  with 
dissimilar  components  (  compare  Chapter  III. A. 3  )  . 


a) 

^4  *  •  4  1  .5  , 

t  *  6 

,  n  =  2 

b) 

Shorthand  notation: 

EXP  (.4)  ♦  EX P  (.  5) 

c) 

Enter  Comment 

Press 

Display 

Initia lize 

z  * 

3 

2 

n 

A 

0 

6 

t 

3 

6 

.4 

.  4 

.5 

/t2 

.  5 

F(t) 

B* 

.  254441493 

calculation  takes  13  seconds 


47 


3.  Find  the  survival  probability  of  a  star.dby-system  wi^h 


one  active  component  and  four  similar  spares. 

a)  K  *  •  3  ,  t  =  7  ,  n  =  5 

b)  Shorthand  notation: 

EX?  (.3)  ♦  EX  P  (.  3)  ♦  EXP  ( .  3)  ♦  EXP(.3)  +  EXP(.3) 


Enter 

Comment 

Pr?  ss 

Display 

Initia  lize 

z 1 

0 

5 

n 

A 

0 

7 

t 

B 

7 

.3 

i 

.  3 

F(t> 

A  • 

. S378738848 

calcula 

tion 

takes  9  seconds 
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PROGRAM  2  :  Procedure 


CASS  I  ;  To  find  the  convolution  of  up  to  four  exponential 
random  variables. 

1.  Use  any  library  nodule. 

se-Partition  (  enter  2  in  the  display,  press  2nd  Op  17  ) 
Read  in  all  four  sides  of  the  magnetic  card. 

2.  Press  2nd  C’  to  initialize. 

3.  Enter  n,  the  number  of  exponential  terms  to  be  convolved 
in  the  display  and  press  A. 

4.  Enter  time  t  and  press  B. 

5.  Enter  hi  and  press  C  (  n  entries  )  . 

REMARK:  Failure  rates,  which  appear  only  once  in  the 
expression,  have  to  be  entered  before  failure 
rates,  that  appear  several  times. 

6.  To  find  the  survival  probability  of  the  system  press  S. 
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PROGRAM  2,  CASE  I  :  Sample  Problems 


(1)  Shorthand  notation 

EXP  (  A/j  )  ♦  EXP  (/it)  +  EXP  (  ) 

Sample  values  :  A4  -  .3  ,  ht  =  .4  ,  t  =  7 
Procedure  : 

Enter  Comment  Press  Display 


Initia lize 


0 


3  n  A 


0 


7  t  B 


7 


.3 

.4 

.4 


hi 

hi 


.  3 
.4 
.  4 


P(t>  E 


. 5363473866 


calculation  takes  14  seconds 
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(2)  Shorthand  notation 


expU„  )  ♦ 

expUj  )  + 

EXP(  /l2  ) 

+  EXP  (  ) 

Sample  values  :  A*  =  . 

2  ,  A  *  = 

.4  ,  t  =  3 

Proc  edure 

• 

• 

Enter 

Comment 

Press 

Display 

Initia lize 

C ' 

0 

4 

n 

\ 

0 

3 

t 

B 

3 

.2 

A, 

3 

.  2 

.  4 

Ai 

2 

.  4 

.4 

Ai 

r» 

.  4 

.4 

hi 

3 

.  4 

F(t> 

E 

.9809746 

calculation  rakes  20  seconds 


51 


(3)  Shorthand  notation 


EXP  (Ka  )  ♦ 

BXP(  /L,  )  ♦  EXP  (  /ii  ) 

+  EX?  (  yit) 

Sample  values  :  yl*  =  .4  ,  hz  = 

.3  ,  t  =  5 

Procedure  : 

Enter 

Comment  Press 

Display 

Initialize  C  * 

0 

u 

n  ft 

0 

5 

t  B 

5 

.4 

/L  c 

.  4 

.4 

Ka  2 

.  4 

.3 

A  l  2 

.3 

.3 

k,\  2 

.  3 

F(t)  E 

. 9029040721 

calculation  rakes  20  seconds 
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(4)  Shorthand  notation 


EXP  (/{,„)  +  EXP  (i2)  +  EXP(/i3)  +  EXP  (  ) 


Sample  values  : 

A/  ^  —  .  1 

!  = 

•  3  f  /£  j  —  .  5  f 

t  =  10 

Procedure  : 

Enter 

Com  ment 

Press 

Display 

Ini 

tia liza 

"  i 

0 

4 

n 

A 

0 

10 

t 

3 

10 

.  1 

i. 

3 

.  1 

.3 

Aji, 

3 

.3 

.5 

.  5 

.5 

A/i 

2 

.5 

F(t) 

E 

.7312684703 

calculation  takes  25  seconds 
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PROGRAM  2  :  Procedure 

CASE  II  :  to  solve  problems  of  the  kind,  that  were 
introduced  in  Chapter  III.B.  . 

1.  Derive  the  system’s  shorthand  notation.  Find  either  the 

-  graphical  representation  or 

-  the  MIX-notation  . 

2.  Dse  any  library  module. 

Re-Partition  (  enter  2  in  the  display,  press  2nd  Op  17  ) . 
Read  in  all  four  sides  of  the  magnetic  card. 

3.  Press  2nd  C'  to  initialize. 

4.  Enter  time  t  and  press  B. 

5.  Repeat  the  following  steps  for  each  path  of  the  graphical 
representation,  i.  e.  for  each  convolution  in  the 
MIX-nota-ion. 

a)  Enter  n,  the  number  of  exponential  terms  to  be 
convolved,  in  the  display  and  press  A. 

b)  Enter  hi  and  press  C. 

REMARK:  Failure  rates,  which  appear  only  once  in  the 
expression,  have  to  be  entered  before  failure 
rates,  that  appear  several  times. 

c)  Enter  p- ,  the  weight  in  the  ith  path,  and  press  D. 

d)  To  find  the  part  of  the  system's  survival  probability, 
that  is  contributed  by  the  ith  path,  press  E. 


6.  To  find  the  survival  probability  of  the  system 
press  2nd  E*. 
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PROGRAM  2,  CASE  II  :  Sample  Problems 

1.  Find  the  survival  probability  of  a  parallel  system 

with  dissimilar  failure  rates  (compare  Chapter  III.B.1). 

a)  >{,4  =  .1  t  A* j.  -  ,2  f  t  =  2 

b)  Shorthand  notation 


EXP  {.  1) 


EXP  (.3)  ♦ 


F(t)  =  MIX[  (.2/. 3)  (  EXP(.I)  ♦  EXP  (.3), 
(.  1  /.  3)  (  EXP  (  .  2)  ♦  EXP(.3)  ]. 
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C) 


Procedure  : 


Enter 

Cod  ment 

Prsss 

Display 

Initia  lize 

^  i 

w 

0 

2 

4. 

B 

2 

2 

A 

0 

.  1 

A, 

n 

.  1 

.3 

Aa  *  Az 

2 

.3 

(.2/.  3) 

** 

D 

.6666666667 

E 

.635793541 

2 

n2 

A 

a 

.2 

At 

r% 

W 

.2 

.3 

At i  ♦ 

r* 

.3 

(.  1/.  3} 

Pz 

D 

.3333333333 

E 

. 304445622 

F(t) 

E* 

.940239163 
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2.  Find  the  survival  probability  of  a  series  system  with 


one  spare  as  introduced  in  Chapter  III.B.2  . 

a)  ^  3  •  3  f  3  •  5  i  t  *  7 

b)  Shorthand  notation 


ZERO 


EXP  (.8) 


F(t)  3  MIX[  (.3/. 8)  (  EXP  (  .  8)  , 

(.5/. 8)  (  EXP ( .  8)  +  EXP  (.8)  ]. 


C) 


Procedure  : 


Enter 

Coo  oent 

Prs  ss 

Display 

Initialize 

Z  1 

0 

7 

4- 

B 

7 

1 

A 

0 

.8 

n 

.3 

(.  3/.  8) 

p« 

D 

.375 

E 

.00  13866  989 

2 

nz 

A 

0 

.8 

/i-t  + 

.8 

.8 

/{.<,  +  ^z 

C 

.8  . 

(.5/.  8) 

Pz 

D 

.625 

E 

.0152536878 

F(t} 

E» 

.0166403867 
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3.  Find  the  survival  probability  of  a  Two-ou+-of-Thre® 


System  as  introduced  in  Chapter  III.B.3  . 

a)  bA  -  .2  /  'll  =  .  4  ,  si 3  =  .  3  ,  t  =  9 

b)  Shorthand  notation 

EXP  (.  9) 


F(t)  *  MIX[  (.2/1.1)  (  EXP  (.9)  *■  EXP  (1.1)), 
(.4/1.1)  (  EXP  ( .  7)  ♦  EXP  (1.1))  , 
(.5/1.1)  (  EXP  (.6)  ♦  EXP  ( 1 .  1 )  )  ] . 
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C) 

Procedure  : 
Enter 

9 

2 

1,1 

.9 

(.2/1.  • 

2 

1.1 

.7 

(.4/1. 

2 

1.1 

.6 

(.5/1. 


Comment  Press 
Initialize  C* 

t 


/i*  +  /tj  ♦ 

+  A*  j 

) 


n  ■ 


A*  +  +  ^>3 

h*  +  «^3 

<>  pz 

n3 

/t,  ♦  X-2  +  X<3 

t 

D  P3 

F(t) 


B 

A 


D 

E 


D 

E 


D 

E 

5* 


Display 

0 

9 

0 

1.1 

.9 

. 1818181818 
.0002624871 

0 

1.1 

.7 

. 3636363636 
.0018043754 

0 

1.  1 

.6 

.4545454545 

.0044892129 

.0065560755 
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